A level coupled to a ID interacting reservoir : A DMRG study 
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The influence of interactions in a reservoir coupled to a level on the width of the filling as a 
function of the chemical potential and the position of the level is studied. The density matrix renor- 
malization group (DMRG) method is used to calculate the ground state of a finite-size interacting 
reservoir, linked to a single state dot. The influence of the interactions in the lead as well as dot-lead 
interactions is considered. It is found that interactions in the reservoir result in a decrease in the 
resonance width, while the dot-lead interaction has an opposite effect. These effects are explained 
within the random phase approximation as an effective change in the inverse compressibility of the 
reservoir, while the dot-lead interactions renormalize the position of the level. 
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There has been much interest in the conduction 
through a ID interacting system, especially in clarifying 
the behavior of a Luttinger liquid with impurities^. Es- 
sentially, it was shown that any impurity will lead to an 
insulating behavior. The resonance conductance through 
a quantum dot coupled to a pair of Luttinger liquid 
leads (see Fig. QJi) was found to produce infinitely sharp 
Coulomb blockade peaks at zero temperature 2 . Thus, 
no level broadening of the dot states is exhibited in the 
measurement of the conduction through that dot. 

Nevertheless, this does not imply that coupling a dot 
to a Luttinger liquid has no effect on the width of the fill- 
ing as a function of the chemical potential. Consider for 
example the arrangement depicted in Fig. QJ). A dot is 
connected to a Luttinger liquid lead, while its occupation 
is measured by a quantum point contact (QPC). Thus, 
the Luttinger liquid acts as a reservoir for the dot, while 
the QPC is used to probe the dots level broadening. In 
such an arrangement, any additional broadening of the 
levels due to the coupling to the reservoir will be seen in 
the shape of the conductance through the QPC. 

The difference between the two arrangements is that 
while the first case (Fig. QJl) essentially probes the en- 
hancement of the backscattering in the Luttinger liquid 
in the vicinity of the Fermi energy, the second (Fig. 
explores the broadening of the level due to coupling to 
states which may be far from the Fermi energy. There- 
fore, one might expect the broadening of the level mea- 
sured in the second arrangement to approach the con- 
ventional Brief- Wigner form, although some signature of 
the interactions is anticipated. 

In this paper we study the broadening of a level cou- 
pled to a Luttinger liquid reservoir (Fig. HJ)). We 
use the numerical density matrix renormalization group 
(DMRG) 3 in order to calculate the ground state of the 
dot-lead system for spinless interacting electrons. It will 
be shown that the level broadening even for an inter- 
acting lead is well described by the Brict- Wigner form. 
Despite that, the interactions in the lead leave a clear 
signature in the broadening of the level. If there is no 
electrostatic coupling between the dot and lead, the inter- 
actions in the lead result in a reduced width of the level. 
On the other hand, such electrostatic coupling leads to 



an increase in the width as well as to a change in the 
levels' energy. 
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FIG. 1: (a) A level (quantum dot) coupled to two Luttinger 
liquid leads represented by the wide lines, (b) A level coupled 
electrostatically to a QPC through which the conductance in 
measured, and to an interacting reservoir, (c) The diagram- 
matic representation of the RPA approximation of the self 
energy. The line corresponds to the lead Green function, the 
black dot to the hopping into the dot and the wiggly line to 
the interaction. 

The DMRG method is used to calculate the orbital 
population in a system consisting of a one-orbital dot 
coupled to an interacting lead. At the first stage in- 
teractions between the dot and the lead are not taken 
into account. The Hamiltonian describing the system is 
therefore given by: 
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where eo is the dot's orbital energy level, V it) is the dot- 
lead (lead) hopping matrix element, and / is the nearest- 
neighbor interaction strength in the lead, a) (a) is the 
creation (annihilation) operator of an electron in the dot, 
and Cj (cj) is the creation (annihilation) operator of an 
electron at site j in the lead. 
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The Hamiltonian H was diagonalized using a finite-size 
DMRG calculation^ for a lead of 150 sites, and n(n) 
curves were calculated for several coupling strengths, V, 
and interaction strengths, /. The lead hopping element, 
t, is taken as 1, in order to set the energy scale. 

Typical results for n{fx) are shown in Fig.|2 As can be 
seen, increasing the interaction strength in the lead (/) 
results in a decrease of the level width, while almost no 
shift in the level position occurs. 
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FIG. 2: (a) The dot's population as a function of the chemical 
potential for different interaction strengths. The full lines rep- 
resents results for lead-dot coupling V = 0.15, and the dashed 
lines for V = 0.3. In both cases the interaction strength / 
takes values between and 2, in jumps of 0.5. (b) The same 
plots for 1 = and 2 are drawn again (V = 0.15 in cir- 
cles, V = 0.3 in triangles) together with theoretical fits using 
Eq.lJBJ (lines). Inset: V^tt as a function of / (symbols) as ob- 
tained by fitting the n([i) curves of V — 0.15 to Eq.©. The 
line corresponds to the dependence according to Eq.JSJ. 



In order to estimate the influence of the interactions 
in the lead on the shape of one should start from 

considering the non-interacting case. The coupling of the 
dot state to the continuum (akin to the Fano- Anderson 
model) may be treated using standard Green function 
technique 5 which leads to: 
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where S(e) is the self energy given by: 
- \V k \ 2 
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6fe are the eigenvalues of the lead, Vjt is the coupling be- 
tween the eigenstates in the lead and the state in the dot 
and S — ► 0. 



For the idealized case, the density of states in the lead 
is constant (i.e., e k = k/Lu, where v is the (constant) 
local density of states, and L is the leads length). The 
coupling is V k = \JajLV (a is the nearest neighbor dis- 
tance), and under these conditions 
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resulting in 3S(e) = ixav\V\ 2 = T/2 and 5R£(e) = 0. 
Thus, one obtains the Breit-Wigner formula: 
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For the tight-binding model given in Eq.Q), e/. = 
—2tcos(ka) and V k = ^2a/Lsm(ka)V, resulting in 
3?£(e) = (V 2 /t)^/l-{e/2t) 2 and 5RS(e) = {V/tfe/2. 
We thus find 
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We shall now turn to discuss the role played by inter- 
actions in the lead. It is well known that the excitations 
in the vicinity of the Fermi energy of any ID interacting 
system should be described as a Luttinger liquid. Nev- 
ertheless, the dot occupation ndotifJ-) is determined by 
contributions from all energies, and the region around 
the Fermi energy does not play a unique role. Therefore, 
one could expect a simple perturbation description of the 
interactions in the lead to suffice. Indeed, the effect of 
the e-e interactions in the RPA approximation on the self 
energy (see Fig. QJ;) may be written as&: 



S(e) = X S°(e) 
where S°(e) is the non- interacting self energy and 
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Here we assumed a constant local density of states 
v (for the tight binding lead v = (a7rf) _1 which ig- 
nores local density of states variations), thus one obtains 
9fE(e) = (V 2 /t{l + I/nt))y/l - (e/2i) 2 and 5RE(e) = 
(V 2 /t(l + I/irt))e/2, corresponding to replacing V 2 in 
Eq.(0 by an "effective" coupling V 2 f j = V 2 (l + I/irt). 

Returning to the results obtained by the numerical 
DMRG calculations, the curves of Fig.|3can be now fitted 
to Eq.© with two fitting parameters - eo,e// and K//- 
It is easy to see that the Brief- Wigner form fits quite 
well even in the presence of strong interactions. The ef- 
fect of interactions is limited here to a decrease of V 2 ^, 
(3) i.e., decrease of T, while the level position eo, e // = (a 
remains constant. The values of V 2 ^ extracted from the 



fit are plotted in Fig. |2mset) and compared with the 
RPA predictions of V 2 ff = V 2 (l + I /id). A rather good 
correspondence is observed. 
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In order to consider interactions between an electron 
occupying the dot, and the electrons in the lead, an ad- 
ditional interaction term should be added to the Hamil- 
tonian (Eq. QJ): 



Hdi = /(flS+ocJci. 



(9) 



DMRG calculations were performed for H + Hdi, and 
the corresponding results (Fig. |3J clearly show a 

change in the resonance width, but also a change in the 
level position, which was absent in the previous case. 
Nevertheless, these results can still be fitted to Eq.© 
with the same fitting parameters - £o,e// and v£ft- As 
can be seen, Eq.© describes this system quite well. 
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FIG. 3: (a) Population of the dot and of the first site of the 
lead as a function of the chemical potential, for V = 0.15. 
Dot-lead interaction was included (Eq.©) taking Idi = I. 
The curves shown are for I between and 3, in jumps of 0.25 
(full lines, dot population), and for 1 = (dashed line) and 
1 = 2 (dotted line) for the lead population, (b) The plots 
for 1 = 0, 1.5 and 3 (symbols) together with the best fit to 
Eq.©. 

The resonance center movement, £o,e//, as well as the 
width, V^ff , that were obtained from the fit can be seen 
in Fig. gj For small values of interaction both grow 
linearly. First lets try to explain the shift in the resonance 
center. As noted, almost no shift was seen for Idi = 0. 
In the presence of weak dot-lead interactions, one may 
approximate Hdi ~ IdWia'a, where ri\ represents the 
average occupation of the first site in the lead. As can be 
seen in Fig.Et, ri\ is not very sensitive to the occupation 
of the dot and may be replaced by its typical value. Thus, 
the energy of the orbital in Eq. © may be rewritten as 
e o.eff = £o + n\Idi- Indeed, this formula fits well the 
numerical results for small values of Idi (as can be seen 



in Fig.|3fa)), for n\ — 0.14. This result agrees well with 
the value n\ ~ 0.15 in the region of the resonance, taken 
from the data of Fig.|3fa). 




o o o o 



o o o o 



0.024 



FIG. 4: (a) £o,eff and (b) V^tt as functions of Idi = I for 
V = 0.15 (symbols) and linear fits for the region I < 0.5 
(lines) . 

A more striking feature is the behavior of V^tf-i i- e i 
the width of the resonance. There is a distinct qualita- 
tive change in the width behavior, compared to the case 
without dot-lead interactions. As opposed to the mono- 
tonic decrease of * , which was demonstrated in Fig. |3 
(inset), Fig.EJa (symbols) shows that V^cf increases with 
J, until a maximal value is achieved around 1 = 2. For 
larger values of interaction decrease in the width is ob- 
served. 

This enhancement of V^tf is associated to the interplay 
between the population of the dot level to the depopula- 
tion of the first site in the lead, ignored in our treatment 
of eo,e//- This leads to a reduction in the effect of the 
dot-lead interaction which results in an increase in the 
width as depicted in Fig. For weak interactions the 
enhancement of V^if ^ s linear. 

Thus, although the Luttinger liquid has a vanishing 
local DOS at the end of the lead in the vicinity of the 
Fermi energy, a level coupled to a ID interacting reser- 
voir is broadened, since all the reservoir states take part 
in the broadening mechanism. Nevertheless, as we have 
seen, the interactions in the reservoir influence the width 
of the resonance. One might gain some insight from the 
following consideration: For the non-interacting case (for 
constant density of states in the reservoir) the width 
is equal to T = 2nav\V\ 2 , which may be rewritten as 
r = 2ir(a/ L)\V\ 2 dN/d/i. The thermodynamic inverse 
compressibility dN/d^j, is affected by the interactions^. 
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Lets consider the compressibility dfi/dN. In the lowest 
order approximation 6 dfi/dN = (Lv)^ 1 + e 2 /C, where 
C is the capacitance of the system. For nearest neigh- 
bor interaction e 2 = al and as usual C ~ L. Therefore, 
dN/dfi = Lvj (1 + avT). Inserting this to the expression 
for r, we get a result similar to the RPA approximation 
results in Eq.||SJ. Although capacitance is proportional 
to the length L of the lead, so is the density of states in 
a ID system, and therefore it has an influence even for 
an infinite lead. 

The conductance through the QPC in the geometry 
described in Fig. QJd is directly proportional to the occu- 
pation of the dots' orbital due to the capacitive coupling 
between the charge of the dot and the QPC (it is as- 



sumed that no tunneling occurs between the dot and the 
QPC) 8 . Thus, in principal, n<iot{lJ), may be read off the 
conductance through the QPC and the effect of coupling 
of the dot to the interacting reservoir can be measured. 

In conclusion, interactions in a reservoir coupled to a 
resonant level leave clear fingerprints on the width and 
position of the resonance. The main influence is the de- 
crease in the resonance width due to a change in the in- 
verse compressibility of the reservoir. On the other hand, 
the dot-lead interaction shifts the resonance position and 
may also enhance the width. 
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